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Les Consequences 
du veillissement d’une canalisation 


Le vieillissement d'une canalisation d'eau potable et sa 
dégradation progressive se manifestent au fil du temps par une 
diminution des performances hydrauliques du réseau, mais 
également par des ruptures qui entraînent différents types de 
dommages. Dans la présente étude, nous abordons le vieillisse- 
ment des conduites elles-mêmes sans traiter des branche- 
ments. Le mauvais fonctionnement hydraulique du réseau suite 
au vieillissement se traduit par : 


Une chute de pression: lorsque la section utile de la 
canalisation diminue à cause de l'entartrage ou de pro- 
tubérances dues à la corrosion, 


Des fuites diffuses: diminuant le rendement du réseau 





Des ruptures: dues à l’action combinée de la corrosion 
et des mouvements de sols. 


Ces différentes détériorations engendrent ri 





des pertes d'eau (et par conséquent une augmen- 
tation de la production), un accroissement des 
dépenses d'énergie lié à l'augmentation des 
temps de pompage, et des interventions directes 
sur le réseau. Divers dommages sont également $ 
engendrés, telle la mauvaise qualité de l'eau, les HAS 
fuites diffuses qui déstabilisent la conduite en vr 
erodant le lit de pose, les ruptures qui provoquent 
des inondations, des coupures du trafic sur les 
chaussées, des coupures d'eau, des dommages 
chez les particuliers et des plaintes des abonnés. ge SS 

Une décision de renouvellement doit donc COUPURE D'EAU 
combiner à la fois une analyse technique des dé- 
faillances et une analyse économique rationnelle 
des choix possibles. 
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Plusieurs approches de renouvellement ont été élaborées par 
l'analyse du vieillissement des réseaux d'eau potable (Elnaboulsi et Alex- 
andre,1996). La flexibilité et la structure mathématique du modèle des ris- 
ques proportionnels nous fournissent l'outil le mieux adapté à l'heure ac- 
tuelle pour analyser les défaillances des conduites d'eau potable et évalu- 
er leur probabilité d'occurrence. Ce modèle décrit les différentes étapes 
de détérioration des canalisations et clarifie l'effet de l'âge sur ce proces- 
SUS. 

Soit T une variable aléatoire positive qui représente le temps de dé- 
faillance d'un sujet appartenant à une population homogène donnée. 
Nous écrivons la fonction de répartition : 


F-(t) = prob(T < t) = | FG) du 
0 


Ou ft est la fonction de densité. La fonction de survie, ST (t) est la proba- 
bilité que T soit supérieure ou égal at: 


Sr(t) = prob(T = t) =1-Fr(t) te[0, +o] 


ST (t) représente le probabilité qu'un evenement arrive apres un certain 
temps. Dans notre cas cette évenement est la défaillance de la canalisa- 
tion. 

Cette fonction est décroissante et varie entre : 


Sr(0) = 1 et lim Sp(t) = 0 


La fonction de densité fT peut s'écrire sous deux formes : 


u Fw | o Fr(t + At) — Fr(t) 
Ir = At=0 At 

1— Srt + At) — 1 +Sr(t) 
Ato At 

—Sr(t + At) + Sr(t) o dSr(t) 
Ato At dt 


fr (t) = ar) = lim mee 


_ jj Di ee 
pee eee O 
re At 
i AER MES ar À 
T Atao Ñ 
= jm PRE ST <E FA) +11 
rue At 


ob(t <T<t+At) 
Mr ) 
At—0 At 


La fonction de risque h(t) spécifie le taux de défaillance instantané 
à T = t des sujets étudiés, sachant qu'ils n'ont pas eu de défaillances 
jusqu'à t et en ont connu une entre t et t+At. h(t) s'écrit : 


bt <T<t+At/T>t 


prob(t<T<t+NNTET) 
Ato At x prob(T = t) 
_ probe <T < t + At) 
— prob(T = DA Ato At 


_ fr(t) 
_ Sr(t) 





La fonction de risque peut être formulée de la façon suivante : 


FO dS 1 dinS;(t) 
SO dt “SO dit 





h(t) = 
Si: Sr(0) = 1, Alors, Sp(t) = exp (— f; h(u)du ) 


Et la fonction de densité peut alors être exprimée par la relation : 


fr(t) = h(Dexp Ef, h(u)du ) 


Nous étudions l'effet des variables explicatives sur le temps de dé- 
faillances. 

Ces variables peuvent étre les ruptures antérieures et l'histoire de 
la maintenance, les caractéristiques de la conduite (matériau,...), des vari- 
ables exogenes comme celles de l'environnement (sol, activité,...). Plu- 
sieurs études de régression ont essayé de déterminer la relation entre le 
temps de survie t et un certain nombre de variables Z. Les modeles les 
plus utilises dans la pratique sont le modele de Weibull (paramétrique), le 
modèle exponentiel (paramétrique) et le modèle de Cox 
(semi-paramétrique). Ce dernier émet l'hypothèse selon laquelle la proba- 
bilité conditionnelle de défaillance à n'importe quel moment est le pro- 
duit d'un premier terme dépendant du temps, reflétant le processus de 
vieillissement, et d'un second terme dépendant du « stress » qui s'exerce 
sur la conduite, ce qui s'écrit : 

h(Z,t) = ho(t) x e7? 
avec : 

hO : fonction de risque de base, arbitraire et non spécifié, ne dépend que 
du temps 
Z : vecteur des caractéristiques du sujet. 
b : vecteur de paramètres propres à chaque variable. 

hO(t) représente le processus de vieillissement des canalisations 
de manière indépendante du phénomène de « stress » auquel celles-ci 
sont soumises. Le facteur stress inclus dans eNZb, est supposé avoir un 
effet multiplicateur sur le taux de risque. 





Le modele des risques proportionnels nous permet de calculer la 
probabilité de chaque défaillance successive en fonction de plusieurs 
paramètres : l'âge de la conduite, son diamètre, la corrosivité du sol, la cor- 
rosivité de l'eau, l’année de pose, le nombre de défaillances déjà subies, ... 

À partir de ces probabilités conditionnelles, nous avons élaboré 
une loi de vieillissement. 

En effet, les différentes fonctions fournies précédemment, c'est- 
à-dire S_T(la fonction de survie à l'instant t), h_t (la probabilité instanta- 
née ou la fonction de risque) et f_t (la fonction de densité) nous ont 
permis de calculer le nombre de défaillances probable à un instant déter- 
miné. 

Si nous considérons que la date de pose t_0=0, alors nous pouvons écrire 
la probabilité d'avoir une défaillance : 
P[N(t)=1] = prob [avoir une et une seule défaillance entre 0 et t]. 
= prob [avoir une défaillance au temps t_1 et ne pas en 
avoir entre t1 et tl. 

= prob [(ne pas en avoir entre 0 et t1 )X (en avoir une en 

t1 )X(ne pas en avoir entre t1 et t sachant qu'il y en a eu une en ti )]. 

=prob [ne pas en avoir entre O et t1] x prob [en avoir en t1 
sachant qu'il n'y en a aucune entre 0 et t1 ] x prob [ne pas en avoir entre t 
et t1 sachant qu'il y a eu une 1ère défaillance en t1 ] 


« Ne pas en avoir entre 0 et t1 » correspond a la fonction de survie S1 (t1). 
« En avoir une en t1 sachant qu'il n'y en a eu aucune entre 0 et t1 » est la 
fonction de risque h1 (t1). 
« Ne pas en avoir entre t1 et t sachant qu'il y en eu une en t1 » est la 
fonction de survie 52 (t-t1) 

Mathématiquement, nous Ecrivons cette probabilité, pour t1 varie 
entre 0 ett, d'avoir une défaillance de la façon suivante : 


PINO =1]= > Sht) = ) (Es) S(t) 


£,=0 E, =0 





Calculons maintenant la probabilité d'avoir deux defaillances entre 


Dett: 
Pour bien comprendre d'ou viendra cette probabilité, on peut se baser sur 


cette illustration : 


AA 


0 t1 t2 t 
-Ne pas en avoir -Ne pas en avoir entre -Ne pas en 
entre 0 et t1 t2 ettl avoir entre t et 
-En avoir une en t1 -En avoir une en t2 t2 sachant 
sachant quil nyen sachantqu'ilyenaeu qu'ily en eu 
a eu aucune entre une en t1 deux en t2 et 
O et t1 t1 
“f1 (t1)” “f2 (t1-t2)” “S2 (t-t2)” 


Donc: 


P[N(E) = 2] = > Fults) y fa — t1) Salt —t2) = y y fi (Ca falls — ty) S(t- 13) 


ES E,=0t,=E, 


En fin, La probabilité d'avoir k défaillances est : 


P[N(t) = k] = y y y e y y fit) fate — t1)falta — ta) ... fie te — tr-1)5k+1 (t — tx) 


t150 ta=tq ta=ta ka Epi Ep Tr 


Cette loi de défaillance nous fournit la probabilité de défaillance instan- 
tanée au temps t. Ceci nous permet de calculer le coût de maintien en ser- 
vice probabilisé et actualisé pour tout événement et a tout instant t. 

Comme nous savons le coút de maintien se change en fonction du 
temps, pour évaluer cette changement, on utilise le taux d'actualisation 
ou de capitalisation. Sion note: 

Cm le coút de maintien en service ; 
Ca le coút d'une défaillance ; 
r le taux d'actualisation ; 


avec: q — (Ca 
m (tr) 








cette dernière relation pour une défaillance entre 0 et t. 


Pour k défaillances entre0ett: Cn = yr En avec: tn=t. 


Soit P(C_m (t)) la probabilité d'obtenir le coút de maintien en service de 
i défaillance entre 0 et t. 
On Constate que : P(C_m (t))= P(N(t)) 


nc: 
E(Cn(e)) = EE Cm ltd) P(Cn(t)) = EE POV) = DIL 


k i 
, 1 
_ ca) P(N(E) = Der: Pe 


1 (1+r)"i E 





Application 
A Résultats 


















































Coút de maintien en service Coút de remplacement 
Nom de la rue Longueur en m Coút direct Coúts indirects 
CAE CF Total Unitaire Total 
Avenue des Vosges 1 150 6 000 108 789 25 920 141 844,5 2 000 2 300 000 



































En premier lieu, le modèle de Cox a été utilisé afin d'analyser l'influence des dif- 
ferents facteurs sur le mécanisme de défaillance et de déterminer les différentes 
courbes de probabilité d'apparition des défaillances et de risque. 

Ces calculs ont été effectués à partir du logiciel EGRET élaboré par le SERC (Statis- 
tics and Epidemiology Research Corporation ). Ce logiciel spécialisé dans l'analyse de 
survie et les analyses épidémiologiques, a permis de déterminer les paramètres du 
modèle de Cox (les différentes variables explicatives, la fonction de risque de base....). 

Ensuite, nous avons élaboré un programme en langage Turbo-Pascal pour résou- 
dre l'équations: k i 

E(Cp(t)) =C, 2, P(N(t) = à) Zar 

En fonction de la date de pose et des conditions structurelles, de l'histoire de la 
canalisation, de sa longueur et des facteurs environnementaux, nous avons calculé la 
date optimale de renouvellement en supposant que cette conduite est soumise a des 
conditions défavorables de détérioration : elles sont supposées étre posés avant 1966, 
de diametre 100 mm, dans un sol corrosif et posées dans une rue a trafic élevé. Le coút 
unitaire de remplacement est de 2 000 F le ml, et le coút direct de réparation est de 6 
000 F. Ce méme calcul a été effectué pour différents niveaux de coúts, a l'aide du tab- 
leau précédent, et en prenant comme hypothèse que les canalisations sont soumises 
a des conditions plus favorables d'exploitation (Elnaboulsi, 1993). 

Tout d'abord, nous avons calculé une premiere fois la date optimale de renouvel- 
lement sans coûts indirects, et une deuxième fois en les intégrant. Nous avons mené ce 
calcul pour une durée de 30 ans à partir de la date actuelle. En outre, nous avons testé 
l'incidence d'une augmentation de 1 000 F du cout unitaire de remplacement. 

En fin, on obtient : 


Avec 
coût 


Coût (Francs) 


10 12 14 16 18 EN 22 24 26 28 30 indirect 
Age (années) 
— Coût de remplacement (F) -- Coût de maintien en service (F) ---.- Coût total (F) 


6 8 0 2 #” 1 3.2 2 4 2 2 20 


Age (années) 
— Coût de remplacement (F) -- Coût de maintien en service (F) ----.. Coût total (F) 





avec coüt indirect 


CONCLUSION 


Cette loi de veillisement nous a 
permis de prendre la décision de répa- 
ration ou du changement du réseaux 
d'eau potable au niveau économique et 
technique. 


Cette loi n'est pas restreinte 
pour les canalisations d'eau potable, il 
peut dévelpoper pour les canalisation 
d'assanissement en se basant sur des 
informations commune au canalisation 
deau potables et d’autre seulement 
pour les canalisation d'assanissement. 


Au niveau du Maroc, cette loi peut 
sappliquer pour estimer le temps de 
défaillance,par exemple l'inondation de 
Casablanca l'année dernière, aussi le 
cout de maintenance et de change- 
ment pour tout le temp. 
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